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Abstract. This article studies the smoothness of conformal mappings 
between two Riemannian manifolds whose metric tensors have limited 
regularity. We show that any bi-Lipschitz conformal mapping between 
two manifolds with C r metric tensors (r > 1) is a C r+1 conformal 
diffeomorphism. This result is due to Iwaniec [11], but we give a new 
proof of this fact. The proof is based on n-harmonic coordinates, a 
generalization of the standard harmonic coordinates that is particularly 
suited to studying conformal mappings. We establish the existence of 
a p-harmonic coordinate system for 1 < p < oo on any Riemannian 
manifold. 



1. Introduction 

This article addresses the following question: given a conformal mapping 
between two smooth (= C°°) manifolds having metric tensors of limited 
regularity, how regular is the mapping? This question is very classical if one 
considers isometries instead of conformal mappings [20|, [8] 0] : it is known 
that a distance preserving homeomorphism between two Riemannian man- 
ifolds with C r (r > 0) metric tensors is in fact a C r+1 isometry. We refer 
to [26] for a recent treatment of this topic and related questions based on 
the systematic use of harmonic coordinates. 

There has been considerable interest in the regularity of conformal map- 
pings between domains in Euclidean space. Classical proofs of the Liouville 
theorem show that any C 3 conformal mapping between domains in M n , 
n > 3, is a restriction of a Mobius transformation and therefore C°°. The 
regularity assumption has been reduced in several works. It is known that 
any Wh™ solution of the Cauchy-Riemann system is a Mobius transforma- 
tion (see [12J), and in even dimensions n > 4 this has been improved to the 
optimal case 

oiW ioc 12 solutions p2]. 
For conformal mappings between general Riemannian manifolds, it seems 
that the literature on such regularity properties is less extensive. It was 
proved by Lelong-Ferrand [7] that a C 3 conformal mapping between two 
C°° Riemannian manifolds is C°°. The work of Iwaniec [11] implies the 
natural regularity for conformal mappings: a conformal mapping between 
two Riemannian manifolds with C r (r > 1 not an integer) metric tensors is 
(jr+l regular. s ee a is [Jg] where the regularity of conformal mappings is 
considered under the assumption that the other Riemannian metric is C r , 
r > 0, and the other one is conformally flat. 
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In this paper we give a new proof that any bi-Lipschitz conformal map- 
ping, or more generally any 1-quasiconformal mapping with W, regularity, 
between two Riemannian manifolds with C (r > 1) metric tensors is in fact 
a C r+1 conformal diffeomorphism. Similarly as in [26], this is done via a 
special coordinate system in which tensors on the manifold have maximal 
regularity. However, instead of harmonic coordinates that are useful in the 
regularity analysis of isometries, we employ a system of n-harmonic coor- 
dinates. The existence and regularity of such coordinate systems on any 
Riemannian manifold is the other main result in this paper. 

To explain n-harmonic coordinates, we recall that a function u on a Rie- 
mannian manifold is called p-harmonic (1 < p < oo) if it satisfies the non- 
linear elliptic equation (p-Laplace equation) 

8(\du\ p - 2 du) = 0. 

Here | • | is the norm induced by the Riemannian metric, d is the exterior 
derivative, and 5 is the codifferential (the adjoint of d in the 1? inner product 
on differential forms). A coordinate system is called p-harmonic if each 
coordinate function is p-harmonic. 

The special property of n-harmonic functions is that pullbacks by con- 
formal mappings preserve this class. In this way, any conformal mapping 
can be expressed as the composition of an n-harmonic coordinate chart and 
the inverse of another such chart. Any system of n-harmonic coordinates 
on a manifold with C r metric tensor (r > 1) always has C£ +1 regularity, 
and the regularity of conformal mappings follows immediately from these 
facts. The n-harmonic coordinates seem like a natural tool, and it is an in- 
teresting question whether they can be used for other questions in conformal 
geometry besides the one treated in this paper. 

This article may be considered as a partial complement to [26], and the 
formulation and proof of Theorem 14.41 concerning regularity of conformal 
mappings are very close to the analogous theorem in [26] for the case of 
isometries. The existence of p-harmonic coordinates, or more generally A- 
harmonic coordinates, for any 1 < p < oo is proved in a similar way as 
the existence of harmonic coordinates (the special case p = 2) in [28] by 
perturbing the Cartesian coordinates and solving suitable Dirichlet problems 
in small balls. However, the underlying equation is nonlinear and solutions 
have a priori only C 1+a regularity for some a > 0, so one needs to work with 
the weak form of the equation. A similar argument was given in [23] where 
solutions to a variable coefficient p-Laplace type equation were obtained by 
perturbing special p-harmonic functions. The present work was motivated 
by [E], where regularity of conformal mappings was a key point. 

Notation. The space C r , when r is a nonnegative integer, denotes the space 
of r times continuously differentiable functions. If r > is not an integer, 
so r = I + a for some integer I where < a < 1, then C r is the Holder space 
consisting of all C l functions whose Ith derivatives are a-H61der continuous. 
We will also use the Zygmund spaces C£, see [27]. By W k,p we denote the 
Sobolev space of functions whose weak partial derivatives up to order k are 
in LP. These spaces are well defined also on smooth manifolds. 
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The matrix of a Riemannian metric g in local coordinates is (gjk), its 
inverse matrix is (g^ k ), and its determinant is \g\. We use the Einstein 
summation convention where repeated indices in lower and upper position 
are summed. 
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2. p-HARMONIC COORDINATES 

If (M,g) is a Riemannian manifold, then the p-Laplace operator is given 
in local coordinates by 

5(\dur 2 du) = -\g\- 1 ^d j {\g\ 1 ^g' k {^a a i^'u) s ^a k u). (2.1) 

The main theorem of this section asserts that whenever (M, g) is a Rie- 
mannian manifold with C r metric tensor and 1 < p < oo, then near any 
point of M there exist local coordinates all of whose coordinate functions 
are p-harmonic. 

Theorem 2.1 (p-harmonic coordinates). Let (M, g) be a Riemannian man- 
ifold whose metric is of class C r , r > 1, in a local coordinate chart about 
some point xo € M. Let also 1 < p < oo. There exists a local coordinate 
chart near xq whose coordinate functions are p-harmonic and have C£ +1 
regularity, and given any e > 0, one can arrange so that metric satisfies 
\djk(zo) — 5jk\ < £ for j, k = 1, . . . , n. Moreover, all p-harmonic coordinates 
near xq have C£ +1 regularity. 

Remark. If T^- are the Christoffel symbols of the metric and T k = g^T^, 
it is well known that the coordinate function x k is 2-harmonic if and only if 
T k = 0. More generally, we have that x k is p-harmonic if and only if 

r k = P -^g k *d t logg kk = (2- P ) 9 —i£-. 

This follows by a straighforward calculation (cf. [13]) from the p-harmonic 
equation flZI]) and the identity T k = — |c/ 1 — 1/2 o , i ( 1 1/2 £/ ifc ) - 

Since this theorem is a local statement, it is sufficient to consider a corre- 
sponding result in M n . Locally, p-harmonic functions in (M,g) are solutions 
of the ^.-harmonic equation div^4(rr, Vu) = in a subset of M n where 

A j (x,q) = \g(x)\ 1 l 2 g^ k (x){g ab (x)q a q h ) E ^q k . (2.2) 

We recall a few standard facts on the ^4-harmonic equation. We will later 
consider a family of equations div^4 e (x, Vm) = where A e satisfy uniform 
bounds with respect to < e < 1, and therefore we need to pay attention 
to the constants in norm estimates. The first fact concerns norm estimates 
for solutions of the Dirichlet problem. 
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Proposition 2.2. Let be a bounded open set in M n and let 1 < p < oo. 
Consider a function 4 : ft, x M n — >■ M n such that for some a, /3 > 0: 

x i-> 4(x, £) is measurable for all £ G M n , 
£ I—)- 4(x, £) is continuous for a.e. x G R n , 
X(a?, f ) • £ > a\£\ p for a.e. x G O and for all £ G R n , 
|4(x,£)| < W 1 for a.e. x G O and for all £ G R n , 
(4(x, £) - ,4(a!, 0) • (C - C) > for a.e. x G O and for ^(Gl n ,^ £. 

Given / G W^fi), there is u G W 1 *^) such that div4(x,Vu) = in Q 
with u — f G W 1,p (r2). Moreover, 

IMIw^fn) < Cll/llw^n) 
where C only depends on a, /?, p, and O. 

Proof. The existence of it for a given / follows from a monotonicity argument 
[HI Appendix I]. For norm estimates, note that u — f G Wo' p (fi) may be used 
as a test function which implies 

/ A(x, Vu) • V(« - /) cte = 0. 

Using the assumptions on 4 and Holder's inequality we obtain 

a / |Vu| p dx< / 4(>, Vw) • Vu = / 4(x, Vtt) • V/ dx 
Jn Jn Jn 

< j ftVvf^Vftdx < /SllVtill^njIlV/H^n). 

It follows that 

||V«||iP(n) < — ||V/||i,p(n). 

Also, by the Poincare inequality 

||« - /||£P(n) < C p , n \\Vu - V/||zP ( n) < C\\V f\\ LP(n) . 

The estimate < C||/||iyi,p(n) follows. □ 

The next result concerns higher regularity of 4-harmonic functions, and 
requires further assumptions on 4. We recall that even in Euclidean space, 
solutions of the p-Laplace equation always have C 1+a regularity for some 
a > but are not C 2 in general. However, if one knows a priori that the 
solution has nonvanishing gradient, then the solution is as regular as the 
function 4 naturally allows. We write B r = B(0,r) for the open ball of 
radius r > centered at the origin in W 1 , and remark that the Zygmund 
space CI coincides with the Holder space C r if r is not an integer [27] . 

Proposition 2.3. Let 1 < p < oo and let 4 : Bi X (R n \ {0}) -> W 1 be a 
C 1 function such that for some 6, (3 > 0, one has for x G B\ and £, £ G M n : 

i4(x,oi + i^4 fe (x,oi + iei i^,4 fe (x,oi < /?ier\ (2.3) 
(4(x,o - 4(x,o) • (* - > + icir 2 ic - ci 2 . (2.4) 
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If u G W 1,p (Bi) is a weak solution of div.4(a;, Vu) = in B%, then there 
are constants C > and a G (0, 1) only depending on n, p, 5, (3, and 
such that u is C 1+a regular in any subdomain of B\ and 

If additionally A is C r , r > 1, on B\ x M n \ {0} then u is C£ +1 in any open 
set where Vu is nonvanishing. 

Proof. If p > n, it follows from Sobolev embedding that 

< ll' u llc' 1 -"/P(Bi) — ( -'n,p\\u\\w 1 'P(B 1 )- 

On the other hand, if 1 < p < n we have by |25[ Thms. 1 and 2] and by the 
bounds A(x,£)-Z> S\£\ p , \A(x,Z)\ < /3|C| P_1 that 

IMIl°°(B 1/2 ) < C n>p> s^\\u\\ LP ^ Bl y 

The assumptions on A also imply that for all x G B\ and £ G M n , 

n 

d id A k (x,0hjh k > 2 p - 2 5\£\ p - 2 \h\ 2 for all h G M n . (2.5) 

The C 1+Q estimate is then a consequence of [3 Thm. 2] (see also [19] for a 
stronger result where A is only assumed to be Holder continuous). 

Suppose then that A is C r regular with r > 1, and that ^ in some 
open set {7 CC Si. We first use Lemma fA. 21 and observe that u G W 2 ^ 2 (U). 
Since u is ^4-harmonic in 17, we have in the weak sense 

dj [A j (x,Vu(x))] = 0. 

Since A G C r and Vu G W 1,2 {U), the chain rule for derivatives holds; see 
e.g. [151 Thm. 6.16]. We have 

d Xj A j (x,Vu(x)) + dt. k A j (x, X7u(x))djd k u = 0. (2.6) 

By the assumption that Vu 7^ in U, in any subdomain U\ CC U we have 
|Vu| > e > for some positive e. The inequality (|2.5p yields in U\ 

d^A j {x,Vu{x))hjh k > 2 p ~ 2 5\Vu{x)\ p - 2 \h\ 2 > 2 p - 2 5e p - 2 \h\ 2 . 

This shows that the equation (|2.6p may be interpreted as a linear elliptic 
equation in nondivergence form. For clarity we write (|2.6p as 

a^ k {x)d j d k u = f, (2.7) 

where the coefficients a? k (x) and the function f(x), defined by 

a jk (x) = d^A j (x,Vu(x)) 

f(x) = -d Xj A j (x,Vu(x)), 

are in the Holder class C a with a = min(a, r— 1) G (0, 1). Also, as mentioned 
u G W io ' c (J7). Thus by elliptic regularity u G C 2+a ; see e.g. O Appx. J]. 

If r is large we can bootstrap this argument to show higher regularity. 
Since r — 1 > cr, we can write r = /co + l + <7 + s where ko is a nonnegative 
integer and < s < 1. We will show that u G (7^0+2+0-^ If fc = this 
has already been proved. If ko > 1 note that u G C 2+<T , so a jfc and / have 
C 1+CT regularity and elliptic regularity (Schauder estimates) imply u G C 3+cr . 
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Continuing this, we obtain u G or equivalently u G C r+1 s . But 

a 3k , f G C""- 1 , souG C£ +1 [271 Theorem 14.4.3]. □ 

We are now ready to show the existence and regularity of „4-harmonic 
coordinates. 

Theorem 2.4 (^4-harmonic coordinates). Let f2 be an open set in M. n and 
let 1 < p < oo. Consider a C 1 function i:!ix(l"x {0}) M n that 
satisfies for x G 17 and £, C G M n , 

) = t p_1 .4(x, for i > 0, (2.8) 

(A(x,Z) - A(x,0) ■ (£ - C) > + Kl) p - 2 |e - CI 2 - (2-9) 

Given any point xo G Q, there is C 1 (in fact C 1+a for some a > 0) diffeo- 
morphism U from some neighborhood of xo onto an open set in W 1 such 
that all coordinate functions of U are ^.-harmonic. 

Moreover, if A G C r with r > 1, then U is C^" +1 and any C 1 diffeomor- 
phism whose coordinate functions are ^.-harmonic has this regularity. Given 
any invertible matrix S and any e > 0, there is U defined near xq such that 

\\DU(x ) - S\\ < e. 

Proof. Assume first that S is the identity matrix. We normalize matters so 
that xq = and B\ C f2, and consider the equation 

6i\A{x, Vu) = in Si. 

The ^4-harmonic coordinates will be obtained as the map U = (u , . . . , n n ), 
where each u 3 = v? e solves for e small the Dirichlet problem 

div^4(x, Vn J ) = in B £ , u 3 \QB e = x 3 . 

The assumptions on A imply that the conditions in Propositions 12.21 and 
12.31 are fulfilled. Then for < e < 1, the Dirichlet problem above has a 
solution u 3 G W 1,P (B £ ) and additionally u 3 has C 1+a regularity. We will 
show that when e is sufficiently small, the Jacobian matrix DU = (djU k )™ k=l 
is invertible at the origin. This gives the existence of the required C 1 (in 
fact C 1+a ) diffeomorphism by the inverse function theorem. The additional 
Cl +1 regularity when A G C r for any system of ^4-harmonic coordinates 
then follows from the previous proposition. 

We will write u = u 1 below. Define the dilated coordinates x = x/e, 
and let u(x) = e~ 1 u(ex) and similarly for (p where tp G C^°(B £ ) is a test 
function. We have Vu(x) = (Vu)(x/e), and 

/ A(x, Vu) ■ Vipdx = / A £ (x, Vtt(x)) • V(p{x)dx = 

J B E JB 1 

where 

A 3 £ (x, q) = A 3 (ex, q). 

Note that A £ satisfies the assumptions of this theorem, and consequently of 
Propositions 12.2 1 and 12.3} uniformly with respect to < s < 1. For instance 

(Mx, o - a £ (x, o) . (e - o > (if i + icir 2 ie - ci 2 (2.10) 
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where the implied constants are independent of e. The function u solves 

divA £ (x,Vu) = in B\, u\qb 1 = x . 

In fact we can take u to be a solution of this equation provided by Proposition 
12.21 satisfying 

INIwLpCBi) < h 
\\u\\ c i +01 (b 1/4 ) ~ 1- 

The original function u may be obtained from u by scaling. 
Write u = uq + u\ where Uq = x . We will prove that 

f \Vin\ p dx < £ min{p ' p ' } (2.11) 

JB! 

where l/p + 1/p' = 1. Define 

1=1 (|Vn| + |Vu [) p ~ 2 |V-Si| 2 (S. 
JBi 

The claim (|2.11j) will follow from the estimate 

/<e||Wi|| LP . (2.12) 

In fact, if p > 2 then 

f \Vu 1 \ p dx= [ |Vni| p - 2 |Vni| 2 dx < f (|V«| + |Vu |) p ~ 2 |Vui[ 2 dx = I 

JBx JB X JBi 

and ()2.12p implies (|2.1ip . In the case 1 < p < 2, the Holder inequality 
implies 

( \Vu 1 \ p dx = ( (|Vn| + |Vn |)" p(p " 2)/2 (|Vn| + |Vn |) p(p ~ 2)/2 |Vni| p dx 
Jb l JBi 

2-p 

< yj (|V«| + 2 I p/2 < I p/2 



since ||'"|| vv^pfSx) ^5 !• Again (|2.12p implies (|2.1ip . 
To show (|2.12p we use (|2.10p and compute 

I=f (|Vfi| + |Vu |) p ~ 2 |Vui| 2 <ii 

< / (^foVfi) -A E (x,Vu )) ■ (Vfi- Vtt )dx. 

Since u solves the _4. £ -harmonic equation and uq solves that .Ao-harmonic 
equation (since Vuq is constant), and since u — uq G wj' p (5i) can be used 
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test function, we obtain 

/ (Ae(x, Vu) - A £ (x, Vuo)) • (Vtt - Vtt ) dx 
JBi 

= — Ae(x, Vtto) • (Vu — V«o) dx 

J B 1 

(A £ (x, Vtt ) - Aq{x, Vu )) • (Vtt - Vu ) dx 

< \\A E {-, V«o) --4o(-,Vfi )||i / p'(B 1 )||V{ti|| XP ( Bl ). 

Here, since .A is Lipschitz continuous and since V^o = e\ is constant, we 
have 

\A £ (S:,Vuo) - A {x,Vu )\ = \A{£x,Vu ) - A(0,Vu )\ < e 
in i?i and therefore 

|| A( • , Vu ) - A( • , Vn )|| lp'(Bi) ~ e - 

The estimate (pTTTj) follows. 

Collecting the results so far, we know that 

||VW - Vu \\ LP{Bl) < £ 'v-^ 

and 

||V«- Vn || C a ( B 1/4 ) £ 1- 
We "interpolate" these estimates by using Lemma lA. 11 which shows that 

1 1 Vtt - Vu \\ L °°(B 1/s ) = o(l) ase->0. 

Repeating this for all components of U = (u , . . . , u n ) implies that 

\\DU(0) - Id\\ = o(l) ase^O. 

Thus in particular DU(0) is invertible for small e, so also DU(0) = DU(0) 
is invertible (and arbitrarily close to the identity matrix) for small e. 

The inverse function theorem shows that U is a C 1 diffeomorphism near 
0, and also U is C£ +1 . To obtain this regularity also for the inverse, write 

DiU- 1 ) = (DUoU- 1 )- 1 . 

Since DU G C£ and U~ l G C r (by the standard inverse function theorem 
and the fact that C£ +1 C C r ), and since C£ is an algebra with respect to 
the usual multiplication of functions, this shows that U^ 1 G C£ +1 . 

The proof of the existence of ^4-harmonic coordinates is completed in the 
case where S = I. If S is any invertible matrix, we can choose u 3 to be 
solutions of 

divA(x, Vu j ) = in B £ , u j \ dBe = {Sx) j . 

Repeating the argument results in ^.-harmonic coordinates for which DU(0) 
is arbitrarily close to the matrix S. □ 

The existence of p-harmonic coordinates is an immediate consequence. 
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Proof of Theorem \2.1\ As discussed above, in any local coordinate system 
where the metric has C r regularity, the p-harmonic equation is of the form 
div A(x,Vu) = where A is given by (pT2j) . satisfies (USD and ([2^j) . and 
has C r regularity. The existence and regularity of p-harmonic coordinates 
is a consequence of Theorem 12.41 

Suppose that the metric in original coordinates at xq is represented by the 
matrix Go, and let U be the corresponding „4-harmonic coordinates. The 
metric in p-harmonic coordinates at xq is given by 

DU^^GqDU (0) -1 . 

1/2 

Choosing U so that DU(0) is very close to G , we can arrange that the 
metric in p-harmonic coordinates is arbitrarily close to identity at xq. □ 

Remark. It is also possible to consider .4-harmonic coordinates on any 
Riemannian manifold (M, g). In this case A is a bundle map T*M — > T*M 
mapping each fiber T* Q M to itself, and the .4-harmonic equation is defined 
naturally as 

5A(du) = 0. 
The conditions f|2.8[) and (|2.9p are replaced by 

A(tO = t^AiO for t > 0, 

(AO - ao, e - o > m\ + icir 2 ie - ci 2 , 

where £, rj £ T*M and ( • , • ) and | • | are the Riemannian inner product 
and norm. The existence of ^4-harmonic coordinates in this setting follows 
by a straightforward modification of the previous arguments. We omit the 
details. 

Standard harmonic coordinates have the property that the Riemannian 
metric has maximal regularity in these coordinates; there is no coordinate 
system where the coordinate representation of the Riemannian metric is 
more regular than in harmonic coordinates. The same is true for ^4-harmonic 
coordinates. More generally, we have the analogue of |13l Corollary 1.4] with 
exactly the same proof. 

Proposition 2.5. Let a Riemannian metric g be C r regular, r > 0, in some 
local coordinates (p near a point xq € M. If a tensor field T is of class C s , 
s > r, in the ip coordinates, then it is of class at least CI in any ^.-harmonic 
coordinates near xq. 

Proof. The chart (p is C°° since M is a C°° manifold, and by Theorem 12.41 
any ^.-harmonic coordinates tp near xq are necessarily C£ +1 regular. The 
coordinate representations T v and of T in the two coordinate systems are 
related by the pullback = (y? -1 o ifi)*T,p. Since the pullback of a tensor 
field involves taking first derivatives of the coordinate functions, T is CI in 
the V'-coordinates. □ 

In the next section we study the regularity of conformal mappings by using 
n-harmonic coordinates, that is, p-harmonic coordinates with p = dim(M). 
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Recall that a Riemannian manifold is said to be locally conformally flat if 
near any point there is a smooth coordinate chart where 

9jk = cdjk 

for some positive function c. Any two-dimensional manifold is locally con- 
formally flat due to the existence of isothermal coordinates. In dimensions 
n > 3, a necessary and sufficient condition for sufficiently regular Riemann- 
ian metrics to be locally conformally flat is that the Cotton tensor vanishes 
for n = 3 and that the Weyl tensor vanishes for n > 4 [2j Theorem 4.24]. 

By the following result, coordinate charts which satisfy gjk = c5jk are 
necessarily n-harmonic. 

Proposition 2.6. If (M, g) is a Riemannian manifold, any coordinate chart 
in which Oj& = cSjk is necessarily n-harmonic. 

Proof. If (x , . . . , x n ) is such a coordinate chart, it follows by a direct com- 
putation that 5{\du\ n ~ 2 du) = for u = x l with any I. (It would enough to 
assume that c is bounded and measurable.) □ 

3. Regularity of C 1 conformal mappings 

The first application of n-harmonic coordinates is to the regularity of 
C 1 conformal mappings between manifolds with C r (r > 1) metric tensors. 
This follows from Theorem 12.11 and the well-known fact that the pullback of 
an n-harmonic function by a conformal mapping remains n-harmonic. For 
C 1 mappings this fact is an easy consequence of the chain rule. In the next 
section we will consider the more general case of 1-quasiregular maps that 

1 n 

have only W^ c regularity. 

Theorem 3.1. Let (U,g) and (V,h) be Riemannian manifolds such that 
g, h G C r , r > 1. If : J7 — > V is a C 1 diffeomorphism that is conformal in 
the sense that <j)*h = eg for some continuous positive function c, then is a 
Qr+l diffeomorphism from U onto V. 

Proof. Choose a point xq € U, and let v be n-harmonic coordinates near 
4>(xq) in (V,h) provided by Theorem 12.11 The map v has C£ +1 regularity. 
Define u = 4>*v, so that u is a C 1 diffeomorphism from a neighborhood of 
x Q into W 1 . 

We claim that all components u l of u are n-harmonic with respect to the 
metric g. To see this, note that when c is a positive function and a; is a 
1-form and all quantities are continuous, we have in the weak sense 

Using the chain rule du l = d(4>*v l ) = 4>*dv\ we compute 
5 g (\du l \ n g - 2 du l ) = S c ^ h (\cp*dv l \ n c ^ h cp*dv l ) 

= c n ' 2 V h (c" ^ (c(f dv l , 0* dv l ) r h ) ^ <j>* dv l ) 

= c n / 2 ct>*(S h (\dvX~ 2 dv l )) = 0. 

Thus u is an n-harmonic coordinate system near xq, and by Theorem 12. II it 
has Cl +1 regularity. 



n-HARMONIC COORDINATES AND REGULARITY OF CONFORMAL MAPS 11 



If r > 1 is not an integer, we have already proved that (p = v o u is a 
(jr+l ma p near X q. The case where r is an integer is handled as in [H [26] 
and is based on the following identity for Christoffel symbols T l - k of (U, g) 
and Tj k of (V, h) under changes of coordinates: 

v%d h <p m - t^di^djcp 1 = did j( p m . 

In our case cp = v ~ l °u \s C T+1 ~ £ for any e > 0, so is at least C 2 and the 
identity is valid. This shows that the second derivatives of (p are C r_1 , and 
consequently <p is C r+1 . □ 

4. Regularity of 1-quasiregular mappings 

In this section we consider the case of 1-quasiregular mappings between 
manifolds with C r metric tensors. The main result, Theorem l4.4[ shows that 
such a mapping is always a C r+1 local conformal diffeomorphism. In the 
formulation and in the proof of the theorem we use the theory of quasiregular 
mappings on R n and on Riemannian manifolds. There are several ways to 
define quasiregular mappings on manifolds [101 [12] , and we will follow the 
approach discussed in [16j . Let us first recall the definition of quasiregular 
mappings on M n . 

Definition 1. A mapping <p : £1 — > W 1 is said to be K -quasiregular if it is of 
Sobolev class W^(Q, W 1 ), its Jacob ian determinant has constant sign a.e. 
on O and if it satisfies the distortion inequality 

\W\ftp < * J* ™. 
If (j> is in addition a homeomorphism, it is called K -quasiconformal. 

In the definition Dip G Lf oc (Q,M nxn ), called the weak differential, is the 
Jacobian matrix of ip consisting of the weak derivatives of the component 
functions (p 1 . The norm || • || op is the operator norm and = det (Dip) € 
Ll (O) is the Jacobian determinant of ip. 

It can be seen from the definition that a 1-quasiregular mapping <p satisfies 

Dcp T D(p = cl nxn a.e. 

Here c is some a.e. positive function. Thus 1-quasiregular mappings are 
conformal mappings in this sense. For the theory of quasiregular and qua- 
siconformal mappings on M n we refer to |12[ [21] [29[ [23] 

Since we are considering mappings between Riemannian manifolds, it 
is natural to consider the Riemannian analogue of quasiregular mappings. 
For our purposes the natural class of mappings is the class of Riemannian 
quasiregular mappings introduced in [IB] . 

Definition 2. Let <p : (M,g) -)• (N, h) be a localizable W^(M, N) map- 
ping between Riemannian manifolds with continuous Riemannian metrics. 
In this case, the mapping <p is said to be Riemannian K -quasiregular if the 
Jacobian determinant of <p has locally constant sign and if it satisfies the 
distortion inequality 

\\D<p\\ n g <KDet g {D<p) a.e. 
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If the mapping (ft is in addition a homeomorphism, it is called Riemannian 
K -quasiconformal. 

Above D(j) : TM — > TN E L™ oc denotes the weak differential of (ft. The 
(pointwise) norm ||_D0|| g and the determinant Det 5 (D(f)) of D(ft are invariant 
quantities defined by the formulas 

\m\ g = -^tr(g-^*h) 1/2 , 

Det g (D<j>) = det [g^^h) 1 ' 2 . 

A mapping is localizable if for every p S M there is an open set U con- 
taining p and a coordinate neighborhood V of (j)(p) such that (f>(U) CC V. 
The condition of localizability is satisfied for example if (ft is continuous. 
The Jacobian determinant of (ft has locally constant sign if the Jacobian de- 
terminant of every coordinate representation of (ft is either non-negative or 
non-positive a.e. 

Riemannian quasiregular mappings share the same analytic properties as 
quasiregular mappings on M n , but also take the Riemannian geometry of the 
manifolds naturally into account. In the case that a mapping (ft : (M,g) —> 
(N, h) is Riemannian 1-quasiregular, the mapping satisfies a.e. the equation 

(ft*h = eg, 

where the pullback is defined in terms of the weak differential, for some a.e. 
positive function c. Note also that a differentiable mapping satisfying this 
equation is 1-quasiregular. For details of these statements and for further 
details of the definition of Riemannian quasiregular mappings see [16] , 

We record the basic properties of Riemannian quasiregular mappings in 
the following lemma. We omit the proof and refer to [16J. 

Theorem 4.1. Let a mapping (ft : (M,g) — > (N,h) be Riemannian K- 
quasiregular and assume that g and h are continuous. Then the following 
hold: 

(1) The mapping (ft is differentiable a.e., and at the points where the dif- 
ferential exists, it coincides with the weak differential. The mapping 
(ft can be redefined on a set of measure zero to be continuous. 

(2) Let u e Wj£(N). Then <ft*u = u o (ft is of Sobolev class W^(N) 
and (ft*u satisfies a.e. the chain rule of derivation: 

di(u o (ft) = d a u\<f,di(ft a i.e. d(ft*u = (ft*du. 

Moreover, 

\d(u o (ft)\ n g < n n KBet [Deft) (ft*(\du\l) a.e. 

(3) If the mapping (ft is non-constant, the Riemannian Jacobian deter- 
minant 

Det (Deft) = v 7 det {g' l cp*h) G L} 0C (M) 
of (ft is non-vanishing a.e. 
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(4) If (j) is in addition a homeomorphism, and thus i^T-quasiconformal, 
its inverse is X n_1 -quasiconformal. Also, integration by substitution 
is valid: 



f fo<f>T>et{D<l>)diJig= [ fdfi h . 

JM JN 

Here / is any integrable function on N. 



Note that the lemma essentially states that all the standard formulas for 
derivatives and integration of smooth mappings continue to hold also for 
Riemannian quasiregular (or quasiconformal) mappings. 

We will also use the following auxiliary result. 

Lemma 4.2. Let (j) : (M,g) — » (N,h) be a Riemannian 1-quasiregular 
mapping. Let p £ M. Given any e > 0, there is 6 > such that if U and V 
are coordinate neighborhoods of p and <p{p) with (j>(U) C V and 

\\9jkip) ~ InxnWop < S, \\h jk (4>(p)) - InxnWop < 

then there are subdomains U' CC U and V CC V such that the coordinate 
representation of <f> : U' — ^ V is an (1 + e)-quasiregular mapping on W\ 

Proof. We will use local coordinate expressions for g, h, and 4>. If 5 is chosen 
sufficiently small, there are subdomains U' and V such that 

lr|y s < 1 + £ » 1/2 " d W) £ < 1 + £ » 1/2 <"» 

on U' and V' respectively. Using the fact that 

<f>*h = eg a.e. (4.2) 

we have a.e. 

We also have det (D</>) = c n ' 2 det (h- 1 ^)) 1 ^ det(g) 1/2 a.e. yielding 

\m\ n op < c^\\h-H<p)\\ n i 2 \\g\\ n l 2 <1 1 cac 

h ~ c n / 2 det(/i" 1 (0)) 1/2 det(<7) 1/2 " 

□ 

We continue with the following proposition, which states that the com- 
position of an n-harmonic function and a conformal mapping gives an n- 
harmonic function. The proof is a straightforward calculation that uses the 
weak formulation of n-harmonicity, but the proof is somewhat lengthy due 
to the weak regularity assumptions. 

Proposition 4.3. Let the mapping <p : (M,g) — > (N,h) be Riemannian 
1-quasiconformal and assume that the Riemannian metrics g and h are of 
class C r , r > 1. In this case, the pullback <p*v = v o (f> of any n-harmonic 
function v on TV" is n-harmonic. 
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Proof. Let v be an n-harmonic function on N. We write for simplicity 
u = v o cf) and show that u is n-harmonic. The function v is at least C 1 by 
Proposition 12.31 and thus u € W^'™(M) by Lemma 14.11 The function u is 
n-harmonic if the functional F, 

F(w):= I g{\du\ n ~ 2 du,dw)dV g , (4.3) 
Jm 

vanishes on the space C^°(M) of smooth functions with compact support in 
M. 

It is sufficient to show that F vanishes for all functions of the form / o 0, 
/ S C£°(N). To see this, let w G C£°(M), if = suppui, and assume that 
F vanishes for all functions of the form / o (f), f £ C^°(N). The function 
k := wo(f)'~ 1 is of Sobolev class W ' n (N) by Lemma [4.1l and by the fact that 
k is supported in the compact set 4>{K). 

Let us approximate k by C^°(N) functions k{ — > k on W 1,n (N). By 
Holder's inequality and the fact that by our assumption F(ki o 0) = 0, we 
have 

\F(w)\ < *) 0)1 + 1^0 0)1 

< ii^iIl" ( m) (/ Wfa - fc ) ° 

Moreover, the following calculation is valid by results of Lemma 14.11 
\d((h - k) o 0)\" g dV g = [ Det {Defy dV g 

< n"/ 2 / - fc)|JJ) o 0Det p0) dVg 

= n«/ 2 / Idih-k^dV^n^Wh-kW^,^. 

J<f>(M) V ; 

Since we have \\ki — &||^/i,n(7v) — >■ 0, i — >• oo, it follows that F(w) = 0. 

We use results of Lemma 14 . 1 1 again to evaluate (|4.3[) for w = f o <f>, where 
/ € C c °°(iV). We have 

"-2..,.. „ r . ..M..7T- _ / ,..|„-L\,.. ^ . ^Det(Zty) 

1 



5 (|o!n| g ra -^n, d(/ o 0))dVg = / _ g{\du\ n - A du, d(f o 0)) — ^ ^ g 



'M JM 

= / \du\ n g -\-ig(du,d(f o^))^ . dV h . (4.4) 

J</,(M) D et {D<f>) | -i 

By the remarks about Riemannian 1-quasiregular mappings, (j) satisfies 

4>*h = eg a.e. (4-5) 

Therefore the Jacobian determinant Det (Deft) equals c™/ 2 . By Lemma 14.11 
the Jacobian matrix Dcj) of 4> is invertible a.e. and the chain rule for deriva- 
tives is valid a.e. yielding 

g(du, d(f o 0)) = g^diiv o fidjtf o 0) = (d^) (d j( f> b ) (0 a v| 

= c(Dci ) )?(D4> T h\ <p Dci ) ) i i(D ( j ) ) b j (d a v\^d b f\t) 

= C (h ab d a V d b f) O (/> = C (h(dv, df)) O cj). 
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Here we have also used (|4.5|) in the third equality. By a similar calculation 
we get an analogous expression for g(du,du). Thus, the right hand side 
of flOp equals 

c\%- 2)/2+1 - n/2 \dv\l- 2 h(dv,df)dV h = [ \dv\ n h - 2 h(dv,df)dV h = 0. 

<j>{M) J<t>(M) 
We conclude that u is n-harmonic. □ 

Our main theorem is a characterization of conformal mappings between 
Riemannian manifolds. The characterization is slightly different when n = 2 
and when n > 3. We begin with the latter case. 

Theorem 4.4. Let (M, g) and (N, h) be Riemannian manifolds, n > 3, 
with g,h € C r , r > 1. Let : M — > ./V be a non-constant mapping. Then 
the following are equivalent: 

<p is a Riemannian 1-quasiregular mapping, (4-6) 

</> is locally bi-Lipschitz and <j)*h = eg a.e., (4-7) 

is a local C 1 diffeomorphism and 0*/t = eg, (4-8) 

is a local C r+1 diffeomorphism and 0*/t = eg. (4-9) 



Proof. I. We show that (|4.6j) implies (|4.9p . First we establish some topolog- 
ical properties of <j>. Let p € M and let £/ be a neighborhood of p and let 
« = («!,... ,v n ) be C£ +1 regular n-harmonic coordinates around 4>(p) such 
that <p{U) C V. 

Since 5 and h are continuous, the coordinate representation of <fi satis- 
fies the M n -definition of a quasiregular mapping [16J. We denote the co- 
ordinate representation of <j) simply by <p : U — > V, where U and V are 
open in M. n . By Theorem 8.13.1 of [12], there is an e > such that every 
(1 + e)-quasiregular mapping <fi : U — > V is locally injective. According to 
Theorem 12.11 and Lemma 14.21 we can reduce U and V such that </> :U —}V 
is (1 + e)-quasiregular mapping on K n . Therefore, by Theorem 14.11 and by 
the remarks above, we can assume that <j> : U — > V is a homeomorphism. 

Proposition 14.31 applied to 4> '■ (U,g) — > (V,h) shows that the pullback of 
the n-harmonic coordinates v = (vi, . . . ,v n ) on V by <f> gives n-harmonic 
coordinates on U. Thus, we have an n-tuple, u = (u , . . . n n ), of n-harmonic 
functions on U and we can express <f> : U — >• V as 

(f) = v^ 1 o u. 

In particular, the regularity of n-harmonic functions shows that (j) is actually 
C r+1 regular as a composition of C r+l mappings at least as long as r is 
not an integer. The inverse of (f> : (U,g) — > (V,h) is also Riemannian 1- 
quasiconformal and thus we can apply the same argument to have the same 
result for . The case r € N is handled by the exact same argument as in 
the proof of Theorem 13.11 

II. One trivially has (gSJ) => (gSJ) => (@2D- Finally, a local bi- 
Lipschitz mapping M — > N is localizable as a continuous mapping. Also 
its Jacobian determinant has locally constant sign, since the local degree 
of 4> is constant and equals a.e. (cf. [lj Thm. 3.3.4]). Since locally the 



n-HARMONIC COORDINATES AND REGULARITY OF CONFORMAL MAPS 16 



Euclidean norm and the distance metric d are equivalent, the coordinate 
representation of is a locally bi-Lipschitz mapping on W 1 and thus belongs 
to the Sobolev class W l( £ for all p > 1. Thus (|4.7f) implies (|4.6p . □ 

In the previous theorem we considered both topological and regularity 
properties of Riemannian 1-quasiregular mappings simultaneously. We ob- 
served that even without a local injectivity assumption the mapping actually 
is a local diffeomorphism. If we assume that the mapping is a priori topolog- 
ically a homeomorphism, we have the following statement that is also valid 
(and well known) in two dimensions. The proof is analogous to that of the 
previous theorem. 

Theorem 4.5. Let (M, g) and (N, h) be Riemannian manifolds, n > 2, 
with g, h £ C r , r > 1. Let (j) : M — > N be a homeomorphism. The following 
are equivalent: 

<j> is a Riemannian 1-quasiconformal mapping, 
4> is locally bi-Lipschitz and 4>*h = eg a.e., 
(f> is a C 1 diffeomorphism and cj>*h = eg, 
4> is a C T+1 diffeomorphism and (j)*h = eg. 

Remark. In two dimensions we cannot expect that 1-quasiregular mappings 
are locally invertible, as shown by the map z \— > z 2 in the complex plane. 

Let us conclude by a simple consequence of the previous results. We call 
a mapping (M, g) — > (N, h) conformal if it is a homeomorphism and satisfies 
one (and thus all) of the conditions (|4,6p - (|4.9p . 



Corollary 4.6. Let (M,g) and (N,h) be Riemannian manifolds with C r , 
r > 1, Riemannian metrics. Let (<pj) be a sequence of conformal mappings 
(M, g) — > (N, h) converging uniformly on compact sets to a mapping (j) : 
M — > N . Then is a C T+1 conformal mapping. 

Proof. A conformal mapping is Riemannian 1-quasiregular. It is well known 
(see [16]) that the sequence ((f) j) converges to a Riemannian 1-quasiregular 
mapping (j). By the previous results, is a C r+1 conformal mapping. □ 

Appendix A. Auxiliary results 

We first give a simple interpolation lemma. 

Lemma A.l (Interpolation). Let O be a bounded open set in R n , and 
suppose that / G C a (Q) where < a < 1. Let also K C SI be a compact 
set, let 6 = dist(K,R n \ Q), and let M > 0. If 1 < p < oo and if 

n-\-ap 

\\f\\ LP (n) < S^M, 

then 

\\f\\L^K)<C n , p , a M^\\f\\^ y 

Here = SU P { ^\x-y\^ i e n,x#y}. 
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Proof. Assume that / is not constant (the lemma is always true if / is 
constant), and choose a point xo £ K. Write B$ = Bg(xo). For < 5 < 5q, 
we compute 

LP(B S ) > \\f(x )\\ L P(B s ) ~ 11/ - f{xo)\\ L v{B 5 ) 

1/p 



1/p 



= \f( X0 )\\Bs\ 1 / p -(j B \f(x)-f(x )\Pdx 

> Cn, p ^\f(x )\ - [f] Ca{U) Ujxrdx 

> c n ,p5 n/ P\f(x )\ - C n ^ a 5 a+n ' p M. 



It follows that 



\L°°{K) < C niP 5 n ^ P \\f\\LP(n) + C n!P!a S a M. 

It is enough to choose 

\fhp(Q ) N n+ap 
M 



□ 



The next lemma shows that ^4-harmonic functions belong to the Sobolev 

2 2 

space W,g C whenever their gradient is nonvanishing. The argument is based 
on difference quotients and is well known (see e.g. [18, Theorem 2.103]). 
However, since we could not find a reference for the precise result needed 
in this paper, the details are included below. See also [22l Theorem 2.5], 
where it is proven that ^4-harmonic functions for p < 2 belong to W?'J! even 
without the additional assumption that the gradient is nonvanishing. 

Lemma A. 2. Let u £ W,'*! (f2) be an ^4-harmonic function, where A sat- 
isfies the assumptions of Proposition 12.31 Assume also that Vu ^ in fL 
Then u G W^(Q). 

Proof. First recall that u € C 1+a (Cl) for some a > [5J Thm. 2], so the 
condition Vu ^ may be understood pointwise. We adopt the notation 

A h f(x) = f(x + h)-f(x). 

A function / belongs to W^ P (Q) if and only if / G L p oc (Q) and 

\\A h f\\ LP{u) < C\h\ 

for any [/ CC H and for sufficiently small h £ M. n [6l Theorem 5.8.3]. Thus, 
to prove the claim, we need to show that 

HAhVull^^ < C\h\, 

for any U CC £1 and for any \h\ sufficiently small. 

Since u is ^4-harmonic, we have for any test function tp G W 1,2 (Q) that is 
compactly supported in Q and for any h sufficiently small that 

/ [A{x + h, Vu{x + h)) - A(x, Vu(x))] ■ V(p(x) dx = 0. (A.l) 
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Choose V so that U CC V CC fi. We will use the test function 

if = r] 2 A h u, 

where rj £ C£°(V), < 77 < 1, is a cutoff function satisfying r/ = 1 near U. 
We separate the expression in the integrand of (|A.lj) above into two parts: 

A(x + h, Vu(x + h)) - A(x, Vu(x)) = h(x) + I 2 {x), (A. 2) 

h(x) = A{x,Vu(x + h)) - A(x,Vu(x)), 
I 2 (x) = A{x + h, Vu(x + h)) - A(x, Vu(x + h)). 

For the first term we write 

A k (x, Wu{x + h)) - A k {x, Vu(x)) 
f 1 d 

= / -7- A k (x,tVu(x + h) + (1 - t)Vu{x))dt 

Jo &t 
1 

d (] A k (x, A]yu{x))A h dju{x) dt. 

Here we have denoted 

A]yu(x) = tVu(x + h) + (l- t)Vu(x). 

Let us denote by DA(£) the matrix (d^A k (x, £)), j, k = 1, . . . , n, where we 
have suppressed the x variable from our notation temporarily. Taking the 
inner product with the gradient of the test function if = rfA^u yields 

(DA(A{Vu)A h Vu) ■ Vcp 

= r] 2 (DA(Alyu) A h Vu) ■ V(A h u) + 2r]A h u(DA{A t h Vu)A h Vu) ■ Vr? 
= S\ + S 2 - 

Now, since A/j Vit = VA^n, the S\ term satisfies 

Si > ci^lA^Vnl^lAftVul 2 . 

Here we have used the ellipticity condition (|2.5|) . To estimate the S 2 term, 
we use Young's inequality with e > 0, 

\S 2 \ < 2\DA(A t h Vu)\\ri A h Vu\\A h uVr]\ 

< \DA{A\Vu)\ (^-\A h uVr]\ 2 + e\rj A h Vu 

Let us next estimate the inner product of V92 and the I 2 term in (|A.2p , 
The assumption (|2.3p and the mean value theorem yield 

\I 2 {x) -Vtp(x)\ < \A(x + h,Vu(x + h)) - A{x,Vu(x + h))\\V(p{x)\ 

< c 2 \h\\Vu{x + h)\ p - l \Vy{x)\ 

< c 2 \h\\Vu(x + h)\ p - 1 (\rj(x)A h Vu(x)\ + 2\riA h u(x)Vri(x)\) . 

Until now, we have only used the structural assumptions of A. We apply 
next the additional conditions u € C 1+a and Vii 7^ in Q, so |Vuj > 5 > 
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on V. Using these facts, we first see that 
Si > c 3V 2 \A h Vu\ 2 

\S 2 \<c 4 (^\A h uV V \ 2 + e\ V A h Vu\ 2 ^ 

\h ■ V<p\ < c 5 \h\ (\rjA h Vu\ + \A h uVr]\) . 
Integrating and using the Cauchy-Schwarz inequality then yields 

= / [h+h] - V^dx 
Jn 

> Ci||7/A h Vu|||2 (v) - C 2 e\\r ] A h Vu\\ 2 LHv) - C 3 \h\ \ \r]A h Vu\ \ L 2 (V) 
C 

- ^\\A h u\\\ 2{v) - C h \h\\\A h u\\ L 2 {v y 

Choosing e and h so small that C\ — C 2 e — C 3 \h\ > C§ > and using the 
mentioned fact about Sobolev spaces that 

\\A h u\\ LHv) < C 7 \h\, 

for all h small enough, yields 

II^V^||| 2(c/) < Csl/il 2 . 

By the discussion in the beginning of the proof, this shows the claim. □ 
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